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Abstract 

We study hot nuclear matter in a model based on nucleon interactions 
deriving from the exchange of scalar and vector mesons. The main new feature 
of our work is the treatment of the scale breaking of quantum chromo dynamics 
through the introduction of a dilaton field. Although the dilaton effects are 
quite small quantitatively, they affect the high-temperature phase transition 
appreciably. We find that inclusion of the dilaton leads to a metastable high- 
density state at zero pressure, similar to that found by Glendenning who 
considered instead the admixture of higher baryon resonances. 
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1 Background of the problem 



Following upon the early efforts [1], |2| to treat nuclei using relativistic mean-field 
theory with the exchange of scalar and vector mesons, a great deal of work was done 
[3-8] to study the consequences of this approach for hot nuclear matter and for 
neutron stars and supernovae P, [TO] . Within the last few years it has been stressed 



11, [L2| that since these approaches to the many-nucleon strong-interaction problem 



must ultimately base themselves on quantum chromodynamics (QCD) they should 
incorporate such symmetries of QCD as chiral and scale invar iance and the breaking 



of these symmetries there at the quantum level. Both references |TT], treat these 
issues at zero temperature, and it is our purpose here to explore the consequences 
of incorporating broken scale invariance, or the QCD trace anomaly, at nonzero 
temperatures. We are interested in particular in temperatures approaching those of 
the deconfmement phase transition, namely, the range of about 150 to 200 MeV. We 
employ a scalar field, the dilaton, to represent the effects of glueballs and to mimic 



the QCD trace anomaly, following a method originally introduced by Schechter 15 



We do not here consider the link between chiral features and the dilaton considered 



by Mishustin, Bondorf, and Rho |T2" |: Its main effect is to displace the minimum of 
the dilaton field from its value when chiral symmetry is respected and the location 
of that minimum is irrelevant in nuclear matter. 



The phase transition found with the present model is obviously not the confine- 
ment phase transition, since the degrees of freedom are entirely hadronic. Neither 
is it the phase transition of chiral symmetry restoration, since our model has no 
chiral symmetry. There may indeed be three separate phase transitions in the 150- 
200 MeV temperature range, where the third is the transition studied here: a sharp 
reduction in the nucleon's effective mass, accomplished by the self-coupling of the 
concomitant high density of nucleon-antinucleon pairs. A more extensive model, 
with the missing degrees of freedom, will show whether two or more of these tran- 
sitions might coalesce into, say, a single quark-hadron transition.0 

Moreover, since the chiral transition is the only one associated with an (almost) 

2 For a discussion of the relationship between confinement and chiral symmetry phase transitions, 
see gcj. 



2 



exact symmetry, it is the only one which assures us that we must encounter an 
actual singularity as we raise the temperature and/or the density. As we shall see, 
the phase transition studied in this paper occurs (as a function of temperature) only 
for a limited range of densities. Elsewhere the high- and low-temperature regimes 
are continuously connected. 

Our motivation for studying dilaton effects in hot nuclear matter is to explore the 
consequences of the fundamental scaling features of QCD within a mean-field model 
that has had much success in its nuclear applications The criticism has been 
voiced [JTSJ that dilaton effects are easily overestimated in nuclear applications if too 
small a mass is attributed to the glueball in place of values for it in the expected 
range of 1.5 to 2 GeV. We anticipate that this issue is irrelevant here, because our 
use of the dilaton in nuclear matter does not require the fixing of a value for the 
glueball mass. Moreover, quite small changes in energy balance within the regions 
of phase transition can lead to different qualitative behavior. We find signs of that 
here. 



As this work was being finished, we learned of a parallel calculation being car- 



ried out by Papazoglou [19j which reaches similar general conclusions about the 
role of the dilaton in hot nuclear matter (and also considers extensions to nonzero 
strangeness). That work considers a particular choice of dynamics for the inclusion 
of chiral invariance and the nucleon mass is generated by spontaneous symmetry 
breaking. It is interesting to compare those results with ours since nuclear thermo- 
dynamics are very sensitive to any change in the nucleon effective mass. 



2 Formalism 



Formal features for the handling of hot nuclear matter in mean-field theory have 
been presented frequently [3-5, 7-12] and we here sketch them only very briefly in 
order to establish notation and procedures. Our lagrangian is |L2], |15"f 

l/) + l&>(Td ll (T-V[<T, 0] 



C = ip 



h* 1 (dp + ig u ^ 
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-~(>V - d»cu») (d,u u - cW M ) + l -ml 



( 



(i) 



where ^ is the nucleoli Dirac field and lo^ is the field for the vector meson. We use 
the usual fourth-order polynomial form for the potential of the scalar field a in the 
presence of the dilaton 0, 



yields the usual classical minimum for the field O = A/e 1/4 , where e w 2.718. The 
coupling strengths for nucleon-meson couplings are g a and g u for the scalar and 
vector fields respectively, and the masses of these fields are m a and m u . We also 
require as input parameters here the coefficients g^ and g± in the scalar potential. 
The powers of appearing in eqs. ([[]) and (0) are determined, as usual, by the 
dimensionality of the term in question. 

The mean-field equations in infinite nuclear matter are given by 





The dilaton potential 




(3) 





for the scalar field, 




(5) 



for the vector field with only the time component u present, and 



M P , + - g ^ + AB £ 1+log^ 




for the nucleon. Here the nuclear scalar and vector densities for nonzero tempera- 



tures are 





and 




(8) 
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The factor 4 appearing in eqs. fl7|) and (§) is for nucleon degeneracy. The fermion 
and antifermion distributions are 

± _ 1 
"» ~ l + exp[(6 ± (p)T^)/T] • (9) 

Here the nucleon and antinucleon energies are given by 



e± = \/p 2 + M* 2 ± 9uj uj (10) 

and the nucleon effective mass is 

M*=(J^M + g a a. (11) 

These equations are now solved for given values of the temperature. Once this is 
accomplished, the system energy density is 

+li m J^-)) alpl, (12) 
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and the pressure is 



p = r*f 7^ -Jrm H + s) - v\ a , ® - c m 



3 7 (2vr) 3 Vp 2 + M*' 2 



Thermodynamic properties of hot nuclear matter are most easily discussed in 
terms of the three intensive variables that characterize the mixed phases near a 
transition, namely, the pressure P, temperature T, and chemical potential fi. This 
immediately recommends the use of the grand potential []2D, ^1] 



Ufa T) = In(Ai, T) = -P . (14) 

Near a phase transition we then have thermodynamic equilibrium for that branch 
having highest pressure (minimal grand potential) at given chemical potential and 
temperature. Thus we show our results as functions of chemical potential rather 
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than baryon density so that transition points may easily be read off the pressure 
curves and applied to the other thermodynamic variables. 

At the temperatures of interest, a pion gas coexists with the nucleon gas we are 
studying. Its mean field is zero, and it affects the nucleon gas through scattering. 
Its treatment is beyond the scope of our mean-field approach. 



3 Results and discussion 



We fix parameters, basing ourselves on previous work ||, such that at nuclear sat- 
uration the density is p = 0.153 fm~ 3 , the binding energy per nucleon is E B = 
— 16.5 MeV, and the compressibility is K = 212 MeV when the dilaton is not 
present and K = 234 MeV when it is. The nucleon's effective mass is M* = 0.57 M. 
The parameters of eqs. (|l|) and @ are then g a = 10.1, g u = 13.3, m a = 492 MeV, 
and m w = 795 MeV for the scalar and vector mesons, and g% = —12.2 fm" 1 and 
#4 = —36.3 for the scalar potential, for the case where the dilaton is absent. When 
it is present the scalar meson coupling becomes g a = 9.6 and the "bag" constant for 
the dilaton field is taken to be B = (200 MeV) 4 , as in O]. Q For nuclear matter 
the value of the dilaton mean field </>o, or alternatively of the glueball mass, is not 
required since we express all quantities in the equations of motion in terms of the 
ratio (p/(po. We note that the scalar potential V of eq. (0) used here suffers from the 
problem of "bifurcation," which is well known in this area j|, [|: In order to obtain 
reasonable values for the nuclear compressibility K we are obliged to use negative, 
and therefore destabilizing, coefficients gs and g±. One must keep in mind that our 
Lagrangian is after all that of an effective field theory. 

Figure 1 shows the well-known ||, ^2| gas-liquid phase transition that occurs 
in these mean-field models at low temperatures. Not surprisingly, there is hardly 
a discernible difference between the cases with and without dilaton for these low 
values of T. 

3 The relative ease with which successful parameters are found in this model comes about because 
we have not constrained it to fit properties of finite nuclei |l3|, . We are satisfied with studying 
dilaton effects on a model which describes well zero-temperature nuclear matter, and thus we have 
not explored the parameter space more fully. 
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We now consider the higher-temperature regime where there appear larger quan- 
titative differences between the two cases. In the theory with dilatons, there is a 
first-order phase boundary in the T—fi plane extending from the T axis at T ~ 
185 MeV to a critical point at (T,fi) ~ (170 MeV, 565 MeV). In the theory 
without dilatons, the intercept is at[] T ~ 189 MeV, with a critical endpoint at 
(T, n) ~ (185 MeV, 300 MeV). Figures 2 and 3 show isotherms of pressure and 
energy density as functions of chemical potential. The phase transition is clearly 
visible as a crossing of pressure curves and as a discontinuity in the energy density. 
To characterize the phases we display the effective nucleon mass in Fig. 4. The 
low (T, fi) regime is a normal phase, by which we mean a state with an effective 
nucleon mass close to that of zero-temperature, ordinary nuclear matter. The phase 
at large T or large fi is an abnormal phase in that the nucleon's effective mass is 
much smaller. For T = 100 and 150 MeV, one sees a smooth crossover between the 
two regimes as \x is varied; a phase transition appears in the isotherms when one 
crosses the phase boundary just described. For sufficiently large T the system is 
abnormal for all \i. 



Since there are small changes in the parameters between the models with and 
without dilatons one should not take the differences in temperature very seriously. 
On the other hand, the abnormal states above the transition region show a rather 
smaller effective mass (Fig. 4) in the presence of the dilaton than they exhibit in 
its absence. This is because when the dilaton is involved it provides an additional 
mechanism to that of the scalar field for reducing the nucleon mass above the tran- 
sition region, as may be seen in eq. (|i~T|). The low effective mass that is found at 
and above the transition region is the characteristic feature of the abnormal branch. 
The occurrence of this abnormal feature is a consequence here || of a feedback loop 
whereby at high temperatures many nucleon- ant inucleon pairs may be produced, 
which in turn enhances the scalar field [see eqs. and (0)], thus further lowering 
the effective mass. In the presence of the dilaton, the <p field also participates in 
this loop through eq. (§). As is well known p3| , there is an alternative mechanism 
which produces an abnormal state at low temperatures purely through the increase 
in baryon density. 



4 Due to numerical difficulties we have not really verified that the phase boundary intersects the 
T axis in the theory without dilatons. 
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The most important consequence of the dilaton is in the metastable extension of 
the abnormal phase to low values of \i. As seen in Fig. 2, both with and without the 
dilaton the metastable state persists all the way down to \i = (whence it continues 
into the /i < region since fi — ► — fi is a symmetry of the theory). In the model with 
the dilaton, however, the metastable state crosses P = which implies that a nugget 
of this phase is mechanically stable. This feature was discovered by Glendenning [H] 
in a model without dilatons which included the contributions of a ladder of baryonic 
resonances to the thermodynamics. We can only speculate that inclusion of both 
the dilaton and the baryonic resonances will serve to drive the zero-pressure state 
towards higher fi and hence closer to the phase transition, thus improving prospects 
for observability. We note, however, that the metastable states in Fig. 2, for both 
models, are separated from the true minimum of the grand potential at the same /i 
(the stable phase) by a very low barrier. This does not make us sanguine regarding 
the lifetime of the state. 



We also present curves for the baryon density (Fig. 5) and the dilaton field 
expectation value (Fig. 6). The latter shows a slow decrease as \x increases and then 
a jump to the abnormal state in which the dilaton field drops by some 10%. 



In sum it appears that the incorporation of scale-invariance breaking in a form 
suggested by QCD can have significant effects in the phase structure of hot nuclear 
matter. We hope that this observation may be of particular use in the study of effects 
of a hot hadronic medium on baryon behavior. This has been explored, for example, 
using skyrmions |2l , though thus far without consideration of dilaton effects. Since 
such effects have been found to be important for one- and two-skyrmion systems 



25| it will be important to take the dilaton into account in both its roles in eventual 



studies of baryons immersed in hot matter. 
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Figure captions 

Figure 1. Pressure P versus chemical potential \x at low temperatures T = 8 MeV 
(right) and 10 MeV (left) in the region of the gas-liquid phase transition: (a) no 
dilaton, (b) with dilaton. Dashed segments represent metastable states. 

Figure 2. Pressure P versus chemical potential ji at the high-temperature phase 
transition: (a) without dilaton, at T = 188 MeV; (b) with dilaton, at T = 175 MeV. 

Figure 3. Energy density £ versus chemical potential ji for values of temperature 
spanning the phase transition region: (a) No dilaton. The isotherms represent, 
from bottom to top, T = 100, 150, 186, 188, and 190 MeV. (b) With dilaton. 
Temperatures, from bottom to top, are T = 100, 150, 172, 175, and 200 MeV. 

Figure 4. Nucleon effective mass M* versus chemical potential /i for values of tem- 
perature spanning the phase transition region: (a) no dilaton, (b) with dilaton. The 
temperature values are the same as in Fig. 3, but in the opposite order: top to 
bottom. 

Figure 5. Baryon density p v versus fx: (a) no dilaton, T = 188 MeV; (b) with 
dilaton, T = 175 MeV. 

Figure 6. Dilaton field 0/0 o for T = 175 MeV. 
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